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Abstract. We introduce nonlocal flows that commute with those 
of the classical Toda hierarchy. We define a logarithm of the dif- 
ference Lax operator and use it to obtain a Lax representation of 
the new flows. 



The present work is based on the joint work pQ of the author with 
B. Dubrovin and Y. Zhang. 

The Toda lattice equation is a nonlinear evolutionary equation in- 
troduced by Toda [2] describing an infinite system of masses on a line 
that interact through an exponential force. In suitable coordinates it 
can be written as a system 

(1) ^u n = e v ^-e Vn 

d 

(2) 777 ^Vi = U n - W n _l, 

at 

with n G Z. 

It was soon realized that this equation is completely integrable, i.e. 
admits infinite conserved quantities, can be solved for rapidly decreas- 
ing boundary conditions through the method of inverse scattering [3] 
and admits explicit quasi-periodic solutions by algebro-geometric meth- 
ods jlj-jHl- It has found important applications in many different fields, 
in particular in the theory of Gromov-Witten invariants of CP 1 , where 
the present extension plays a particular role j7j . 

The Toda lattice equation can be seen as the first element of a whole 
hierarchy of commuting flows, the Toda lattice hierarchy. We can write 
all the flows by using the Lax representation |3J 

< 3 > <wr l iTv (L,+1)+ ' L] 
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where L is the following difference operator 

(4) L = A + u + e v A-\ 

Here we are using a notation (see jH]) with a continuous space variable 
x = ne, where e is the lattice spacing and t = et; i.e. u(x) = u n and 
v(x) = v n for x = en. We denote by A the shift operator Af(x) = f(x+ 
e). Given any difference operator A = J2k€Z a k-^- k we denote by A + , 
A_ the positive and negative parts respectively, i.e. A + = J2 k>0 o-k^ k , 
A = A + + A_. The first flow to is simply given by the Toda equations 
CO)-®; as a further example the second flow is given by 

(5) eu tl (x) = -((u(x + e) + u{x))e v(x+e) - {u{x - e) +u{x))e v{x) ) 

(6) ev h (x) = ]^{e v{x+e) - e v{x ~ e) + u\x) - u 2 (x - e)). 

This hierarchy is Hamiltonian with respect to two different Hamil- 
tonian structures j^j; in our continuous formulation the first is given 
by the following Poisson bracket 

(7) {u{x),v{y)} 1 = -{5(x-y + e) - 5(x - y)) 

e 

and {u(x),u(y)}i = {v(x),v(y)}i = 0. The second Poisson bracket is 

(8) {u(x), u(y)} 2 = ke v ^6(x - y + e) - e v ^8(x - y - e)) 

(9) {u(x), v(y)} 2 = -u{x)(6(x - y + e) - 5(x - y)) 

(10) {v(x), v{y)} 2 = -(S(x -y + e)-8(x-y- e)). 

e 

The equations of motion can be written in the following form 

(11) ^. = { ,/, }l= 1 {,/, } 2 

dt q q + l 

where h q = J h q dx and the Hamiltonians are given as traces of powers 
of the Lax operator L 

(12) K = ^yRes(L^); 

given any difference operator A = J2kez ak ^- k ' ^ s res idue is defined by 
Res A = a . 

The presence of a two Hamiltonian structures permits to obtain all 
the Hamiltonians through the Lenard-Magri ^U] recursion relation (fTTj) 
starting from a Casimir of the first Poisson bracket. In our case if we 
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start from the Casimir h-i = u we obtain all the Hamiltonians de- 
fined above. Moreover this procedure guarantees that all the resulting 
Hamiltonians commute among themselves. 

In order to see why and how this hierarchy of equations could be 
extended, let's consider its dispersionless limit that is obtained putting 
e — > 0. It can be shown [H] that the Lax representation of the disper- 
sionless flows is given by 

(13) wr { J^Ty {a+1)+Xh 

in this case £ is a function of x and of the additional variable p 

(14) C = p + u{x) + e vix) p- 1 
and the bracket is 

9BdC - 9CdB 

^ dp dx ^ dp dx 

for B and C functions of p and x. (£>)+ means that only non-negative 
powers of p, in the power series expansion of £>, are considered. 

The dispersionless Hamiltonians and the Poisson brackets are simply 
obtained from their dispersive counterparts (JT2J), (f7|)- (fTUj) by putting 
e — > 0. In particular the same recursion relation as above (jllj) holds in 
the dispersionless case. 

Considering the genus zero approximation of the topological CP 1 
model, in [TT] it was noted that new flows can be added to the usual 
dispersionless flows given above; their Lax representation is 



BC 2 

(16) - r = {-{C\\ogC-c q )) + ,C} 

ot q q\ 



where c q = J2l=i h c ° = 0- ^ ne logarithm of £ must be understood 
in the following way 

(17) log £ = -v + - log(l + up- 1 + e°p' 2 ) + - log(l + ue~ v p + e" V) 

where the first logarithm on the RHS is seen as an expansion in negative 
powers of p while the second one in positive powers of p. 

These flows can be expressed in Hamiltonian form by 
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where the Poisson bracket is the dispersionless limit of (jZJ) and the 
dispersionless Hamiltonians are given by 

(19) hf* = j^— ^Res^- 1 C^\\og C - c q+l )]. 

These Hamiltonians however satisfy a recursion relation that is different 
from the previous one 

(20) {-, l q _ x } 2 = q{; l q } 1 + 2{; VJl 

We would like to briefly mention that in the dispersionless limit all 
the flows can be introduced in a rather different fashion, by using the re- 
lation between systems of hydrodynamic type and Frobenius manifolds 
[12]. The dispersionless Toda Poisson pencil is associated to a Frobe- 
nius manifold characterized by the free energy F = \w L v + e v ; then 
all the Hamiltonians of the system can be obtained by an expansion of 
the so-called deformed flat coordinates of the Frobenius manifold. The 
relationship between the dispersive and dispersionless versions of the 
Toda hierarchy is just a particular instance of the classification pro- 
gram of bihamiltonian integrable hierarchies proposed in ^Hj , based on 
the reconstruction of the whole dispersive hierarchy starting from its 
dispersionless limit. 

Thus we have seen that in the dispersionless case the Toda hierarchy 
has two perfectly well defined sequences of flows all commuting between 
themselves, denoted by the times t q and t q for q > 0. The classical 
dispersive flows corresponding to the times t q defined above reduce, for 
e — > 0, to the corresponding flows in the dispersionless hierarchy; on the 
other hand in the classical dispersive formulation there is apparently 
no flow reducing for e — > to the dispersionless flows corresponding to 
the times t q . 

In analogy with the Lenard-Magri procedure for the first set of 
Hamiltonians one expects to find the second set of flows of the dis- 
persive hierarchy from another Casimir of the first Poisson bracket by 
applying the recursion relation (J2(J|) (this time with the full dispersive 
brackets). The first Poisson bracket admits in fact a second Casimir: 
h_i = v. However the recursion relation (|2U|) fails to work in this case, 
as one can easily check, the reason being that h-\ is a Casimir of both 
Poisson brackets. This phenomenon is referred to as "resonance" of 
the Poisson pencil. 

We can however introduce an ansatz for the first nontrivial Hamilton- 
ian h ; it was given in [T3] and [T3] and is such that the corresponding 
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flow coincides with the x-translation 

(21) h = uA(A-l)- 1 ev x . 

This Hamiltonian is nonlocal, since it contains the inverse of the dis- 
crete derivative; this can be written as a formal series in e 

(22) (A-l)-V = X;ff (e9x)fci; 

fc>0 

where the Bernoulli numbers B>k are defined by = J2k>o "ff 3 ^- 
Starting from this ansatz one can define all the Hamiltonians h q us- 
ing the recursion relation (}2T)j) and show that they commute between 
themselves and with all classical Toda flows. 

To have an explicit form for these new flows, it is important to 
find their Lax representation. By considering the dispersionless Lax 
representation ()16|) one expects that it will be necessary to introduce a 
logarithm of the Lax operator L; such operator can be defined by the 
dressing formalism. It is well known JE] that one can write L as the 
dressing of the shift operators A and A -1 

(23) L = PAP' 1 = QA^Q- 1 
where 

(24) P = J2PkA- k Po = l 

fc>0 

(25) Q = ^g fc A fc . 

fc>0 

By substituting in the definition the functions pk, qu can be found 
in terms of u, v. 

Noticing that A = e tdx one is led to define two different logarithms 
in the following way 

(26) log + L = PedP- 1 = ed + PeP' 1 

(27) log_ L = -QtdQ- 1 = -ed - QeQ' 1 . 

These logarithms are both differential-difference operators, due to the 
presence of ed; since we want to write an expression like (JT5j) and we 
need to make sense of the (•)+ part, we would like to have a purely 
difference operator for the logarithm, that we define by 

(28) log L = X - log + L + X - log. L = -\(PxP~ l ~ QxQ~ l ); 
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in this definition the derivative drops out and we get a difference op- 
erator of the form 

(29) logL = ^«vV fc . 

fcez 

One would like to express the coefficients Wk in terms of the vari- 
ables u, v. This is indeed possible (see |T] for a proof); essentially, by 
definition, all the previously defined logarithms commute with L, e.g. 

(30) [log + L,L]=0; 
inserting 

(31) log + L = td + 2 ™* A * 

k<-l 

into (j3()|) one can solve it recursively and express the coefficients Wk for 
k < —1 as formal power series in e with coefficients that are differential 
polynomials in u and v. The analogous expression for log_ L gives the 
coefficients Wk for k > 0. First few examples are 

(32) w-i = ^(A-irV, 

(33) w = ^A(A - l)-^ 

(34) Wl = iAe-^A-l)" 1 ^. 

Observe that this is in general not possible for the coefficients Pk, of 
the dressing operators. 

Using this definition of logarithm of L we can give the Lax representa- 
tion for the new flows, which is formally analogous to the dispersionless 
Lax representation 

BT 2 

(35) e- r = [A q ,L] A q = -[L«(logL-c q )) + ; 

here A q is a difference operator of infinite order; equivalently we can 
use the operator 

(36) A q = -\U{\ogL - c q )) + - ^(log_ L - c q ) 

q\ q\ 

that also contains a differential part but is of finite order; it gives the 
same Lax equations since it differs from A q by a part that commutes 
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with L. For example, the first two Lax operators are given by 

(37) A = ed 

(38) A x = A(ed - 1) + A(A - l)~ l eu x + u(ed - 1) 

+e v (ed + 1 - (A - l)~ 1 ev x )A~ 1 . 

As expected the first one corresponds to the x-translations while the 
second gives the first nontrivial extended Toda flow 

(39) eu k = (A-lX-e^A -1 - l) - ^) - 2(A - l)e" 

u 2 

(40) ev h = ((A- 1 -l)- 1 ev a ,)(A- 1 -l)tt + eu 3! (A- 1 «) 

+A~ 1 eu x + eu x + 2(A _1 - l)u. 

We can also give an explicit expression for the Hamiltonians in anal- 
ogy with the dispersionless case 

(41) h q = j-^—Res[L« +1 (logL - c q+1 )]; 

these are exactly the Hamiltonians defined above by the recursion re- 
lation and the ansatz ([21)1 . up to total derivatives. 

The new flows of this extended Toda hierarchy are nonlocal but nev- 
ertheless have many of the nice properties that are expected from a 
completely integrable system, such as the existence of multisolitonic 
solutions; for example, the evolution of a simple soliton under the non- 
local flows ([3~9* J) -PD ]) is given by 

cosh(i((a; + AifOlogZ! + + 1/ Zl ))) 

u{x, tx) = (A — 1) z = 

cosh(-((x + Xxtxjlogzx + tx(-zx + 1/zi)) - logzi) 



f(x,ti) = (1-A 1 ) 2 log[2cosh(-(x + Aiti)logzi + -ti(-zi + l/zi))] 

where Ai = Zx + ^f 1 and is a parameter. As usual this is simply 
obtained by Darboux transformations of a particular constant solution. 

One expects that also the algebro-geometric quasi-periodic solutions 
should have a nice behavior under these flows; this problem is however 
still under consideration. Further generalization for the multicompo- 
nent case is in progress. 

All the proofs of the above statements will appear in pQ. 
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